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ABSTRACT 

Asteroids are frequently colliding with small projectiles. Although each individual 
small collision is not very important, their cumulative effect can substantially change 
topography and also the overall shape of an asteroid. We run simulations of random 
collisions onto a single target asteroid represented by triaxial ellipsoid. We investigated 
asteroids of several hundred meters to about 18 km in diameter for which we assumed 
all material excavated by the collision to escape the asteroid. The cumulative effect 
of these collisions is an increasing elongation of the asteroid figure. However, the 
estimated timescale of this process is much longer than the collisional lifetime of 
asteroids. Therefore, we conclude that small collisions are probably not responsible 
for the overall shape of small asteroids. 

Key words: minor planets, asteroids: general - methods: numerical. 


1 INTRODUCTION 


Collisions between aster oids are a very important process , 
which affects their spins (iHarr is fc Buri3ll979l, Harrisll 19791) . 
size distribution (lDohnan£j|_ 196^ Cel]ino_et_^ 199ll) and 
also t heir shapes ( ICat ullo et al.lll984l . iFarinella fc Zappalal 
I 1997 I . iTanga et al.l 20091 ). Subcatastrophic collisions, how¬ 
ever, are not very well described yet and in evolutionary 
models of asteroid populations, they are only roughly ap¬ 
proximated. 

The role of small collisions on the overall shape 
of asteroids was stud ied in two-dimensional model of 
iRonca fc Furlon'3 ll 19791 ). They used similar assumptions in 
the construction of their model as we do and for the initially 
elliptical asteroid model they ob tained increasing elongation 
of its shape. Later. iHarri 3 (Il99d) mentioned this effect when 
discussing collisional evolution of the spin of nonspherical 
asteroids. 

Numerical investigatio n of asteroid shape evolution by 
collisions was presented by Korycanskv fc Asphaud (1200 ll ) 


and iKorvcanskv fc Asphaug ( 20031 ) (hereafter KA03). Their 
model is based on excavation, orbit and reimpact of the 
ejecta from the cratering event and also on the conserva¬ 
tion of the maximum angle of repose of the material on 
the surface of the model asteroid. KA03 conclude that the 
final shape is mostly oblate ellipsoid indepedently of the ini¬ 
tial conditions. They also calculated the timescale of such 
impact-induced reshaping of the asteroid to be more than 
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an order of magnitude longer than catastrophic disruption 
timescale. 

ISzabo fc KissI (l2008h report about shape distribution of 
asteroid families derived from multi-epoch photometry of 
the Sloan Digital Sky Survey Moving Object Catalog. They 
conclude that younger families seem to contain more elon¬ 
gated asteroids, while older families contain more rounded 
asteroids. In their opinion, this is in agreement with the pre¬ 
dictions of an impact-driven evo lution of the asteroid shapes 
(KAOS). iDomokos et al.l (120091) studied the role of micro¬ 
collisions on the overall shape of asteroids in their averaged 
continuum abrasion model. Their asteroid models evolve to 
prolonged shapes that possess large, flat areas separated by 
edges. 

In the present paper we want to discuss the effect of a 
large number of subcatastrophic collisions on the shape of an 
asteroid. This problem caught our interest when construct¬ 
ing an evolutionary model of asteroid rotations influenced 
by subcatastrophic collisions focused on tumbling asteroids 
dPravec et al.ll2()()^ . We noticed that our model asteroids 
subjected to collisions became ever more elongated than at 
the beginning of the simulation. Since we wanted to ensure 
it was not just a numerical issue of our model, we started to 
investigate this effect. 

We give a detailed description of our model and the 
method used in Section Section [3] contains an estimate 
of timescale of the increasing elongation process. Numeri¬ 
cal results are given in Section [3] In Section [S] we discuss 
some implications and caveats of the present research and 
conclusions are given in Section |6l 
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2 INCREASING ELONGATION OF THE 
MODEL ASTEROID 

An obvious effect of a small collision is a formation of im¬ 
pact crater on the surface of an asteroid. There is, how¬ 
ever, a cumulative effect of a large number of such cratering 
events as well. In our analytical model of subcatastrophic 
collisions between a target asteroid and a large number 
of small projectiles, we noticed that the shape of an ini¬ 
tially elongated asteroid becomes gradually more elongated. 
More exactly, it was changing the shape of a dynamically 
equivalent equal-mass ellipsoid. Its se miaxes are defined as 


equivalent equal-mass ellipsoid, its se n 
llKaasalainenll200ll . iPravec et al.ll2014| j 
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and 


la = h/Is Ib = I 2 / Is ■ 


( 1 ) 


where 7i, I 2 and Is are the principal moments of inertia of 
the target asteroid and la and lb are normalized principal 
moments of inertia of the dynamically equivalent ellipsoid; 
/c = l. 

In our model, the target asteroid was represented by 
a triaxial homogeneous ellipsoid. We ran a series of con¬ 
secutive collisions with the target, each forming an im¬ 
pact crater on the surface of the ellip soid according to the 
scalin g laws of iHolsappl^ (liggsl l and iHolsapple fc HousenI 
l|2007l l. with the angular mom entum transfer effic i ency b ased 
on laboratory experimen t s of Yanagisawa et al.l (Il996l l and 
lYanaeisawa fc Hasegawei (l20od l and updating the inertia 
tensor of the target asteroid after every single collision. 

The impact speed was 5kms“^, which is the median 
impa ct speed in the cen tral regions of the Main Asteroid 
Belt l|Bottke et al.l [200^. The geometry of collisions was 
isotropic with respect to the target surface (the probabil¬ 
ity of a collision from any direction was the same) with one 
exception: We only allowed collisions for which incidence an¬ 
gle (the angle between the impact velocity vector and the 
normal to the surface at the impact site) was less than 75 
degrees. We simply wanted to avoid almost tangent impact 
geometries, since the formation of impact craters for such 
impact geometries is not very well understood yet. 

The target was in a principal axis rotation state before 
the collision. Initial rotation period was chosen randomly 
according to the observed distribution of rot ation periods 
of sm all Main Belt/Mars Crossing asteroids (|Pravec et al.l 
I2OO8I . updated on 2014-04-20). Size of the p rojectile was 
rand om according to a power-law distribution JBottke et al.l 
l2005h . 

After every collision, rotation state and shape of the tar¬ 
get asteroid changed. However, in this paper we only focus 
on a change of the shape and we quantify this change as fol¬ 
lows. The collision formed an impact crater on the surface of 
the ellipsoid. We therefore calculated an inertia tensor of the 
mass that fills the crater and subtracted it from the inertia 
tensor of the ellipsoid. Then, we updated lengths of axes of 
the ellipsoid using the values of the dynamically equivalent 
ellipsoid. We also reset the shape of the target asteroid to 
the triaxial ellipsoid (erased the impact crater) keeping its 
density on the initial value (which lead to decreasing size of 
the target asteroid). 

In our simulations, we only investigated the effect of 


subcatastrophic collisions. This corresponds to the calm 
period of an asteroid’s life between larger shattering or 
even dispersal events (shattering causes large scale dam¬ 
age to the asteroid, but most of its pieces remain gravi¬ 
tationally bound; disp ersal event reduces the target to one - 
half its original mass, iMelosh fc Rvanlll997l . iHouse nl l2009ll . 
We used two criteria to recognize a subcatastrophic colli¬ 
sion. First, we compared the specific energy of the colli¬ 
sion wit h the shattering crite r ion. I n our calculations, we 
followed IStewart fc Leinhardtl (|2012| ~) , who derived the dis¬ 
persal criterion in gravity regime of collisions in their nu¬ 
merical simulations. We calculated the sha t tering criterion 
as 1/4 of their value according to [Housed (l2009l l. Second, 
we calculated the maximum size of an impact crater for a 
gi ven size of the target asteroid a ccord ing to the relation 
of lBurchell fc Leliwa-Koovstyhsl^ (120101 1. The collision was 
subcatastrophic if both these conditions were met. 

There is a set of mechanical properties of the target and 
the projectile, which can be varied in a physically plausible 
range. We did th i s to s ome extent in our previous paper, 
iHenvch &: Prav^ (l2013l l. and we found that these values af¬ 
fect the size of an impact crater formed on the surface of 
the target asteroid. Therefore, they affect the timescale of 
increasing elongation of the target’s shape. In the present 
research, we have chosen some plausible values of those pa¬ 
rameters and held it constant throughout our simulations 
as specified in Table [T] The density of the target and of 
projectiles was the same, 2gcm“®. 

Since we are mostly interested in the small asteroid pop¬ 
ulation (bodies of several kilometers in diameter), we ran 
our simulations for asteroids of mean diameters 0.6-17.6 km 
(see Table [ 3 ). The upper limit was determined by the fact 
that almost all material excavated by the hypervelocity im¬ 
pact of a projectile is ejected from the asteroid for asteroid 
diameters of less than about 18-20 km according to scal¬ 
ing laws (iHousen fc Holsappld l201l| j. This was calculated 
by using the scaling l aws f or ejecta speeds as described by 
lliousen fc Holsapplel ll201lll and comparing it with the sur¬ 
face escape speed. 

The mean diameter of the target is 

A=2(a6c)l/^ (2) 

where a, b and c {a ^ b ^ c) are the semiaxes of the tri¬ 
axial ellipsoid representing the target asteroid. In our sim¬ 
ulations, we concentrated on two effectsfhe shape and the 
size dependence of the collisional erosion. Therefore, we ran 
longer series of 1,000 collisions for various target asteroid 
sizes (see Table [3 and fixed initial axial ratios of a/c = 2.0 
and bjc = 1.4. For each target size we ran three separate 
runs to exclude any numerical artifacts. Then, we ran shorter 
series of about 110 consecutive collisions, where we fixed 
c = 0.5 km and varied the initial semiaxes ratios (see Ta¬ 
ble [ 3 . 


3 TIMESCALE ESTIMATES 

The importance of erosional process for evolution of aster¬ 
oid shapes is determined by the timescale of this process. We 
calculated the timescale in two different ways and obtained 
very different results. We also have the timescale indicated 
by the numerical simulations as explained in a following sub- 
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Table 1. Scaling and material constants for non-porous 
target. Ki and K 2 are experiment-based coefficients, 
and u are the scaling parameters, y is a materi al strength 
and K r and are the shape constants, cf. fHolsappld 

1I2003I') . 


Ki 

K2 

h 

u Y (MPa) 

Kr 

Kd 

0.095 

0.257 

0.55 

0.33 1 

1.1 

0.6 


Table 2. Effective diameters of model asteroids in long se¬ 
ries of collisions and limiting projectile sizes. The fixed initial 
semiaxes ratios of the target were ajc = 2.0 and h/c= 1.4, Dt 
is the target mean diameter, Dp^max is the largest projectile, 
^p,min is the smallest one. 


c [km] 

0.2 

0.35 

0.5 

1.5 

2.5 

4.25 

6.25 

Dt [km] 

0.6 

1.0 

1.4 

4.2 

7.0 

12.0 

17.6 

-^^p,max 

5.5 

10.5 

16.0 

78.0 

150.0 

291.0 

437.0 

Dp,min H 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 


section. We discuss the differences between the two timescale 
calculations at the end of this section after their description. 

3.1 Timescale of the erosion from the simulations 

The timescale of the increasing elongation of the asteroid 
figure can be estimated directly from the simulations. We 
filtered size of the projectiles for a specific size of the target 
to keep the collisions subcatastrophic (see Table [2J. How¬ 
ever, we also saved information about discarded projectiles, 
whether too small or too large. Therefore, we have the rela¬ 
tive timescale of the increasing elongation in comparison to 
the catastrophic disruption of the target asteroid caused by 
projectiles larger than upper size limit. 

In all our simulation runs with 1,000 subcatastrophic 
collisions, the target asteroid was hit be the projectile larger 
than upper size limit after a few subcatastrophic collisions. 
Therefore, we estimate the timescale of the process to be 
at least two orders of magnitude longer than the collisional 
lifetime of the target asteroid. 

3.2 Ejected mass measure 

The first and probably more reliable and solid calculation 
compares the erosional effect of projectiles in terms of exca¬ 
vated and ejected mass. As noted above, practically all mass 
excavated by the projectile escapes the target asteroid for 
sizes used in our calculations. The total eroded mass by the 
population of projectiles can be calculated as 

r^p,max 1 

Me. = / -pi{Dt + D^fMAV^n{D^)p^, (3) 

^p,min 

where pi is the intrinsic probability of collision, Dt and Dp 
are the target and the projectile diameters, respectively. At 
is the timescale we are looking for, Vp is the volume of the 
projectile, A is a factor describing the excavated material 
volume, n{Dp) is the distribution function of the projectile 
sizes and Dp,min and Dp,max are the minimum and maximum 
projectile size, respectively. 


We define the timescale of erosion process as a time 
interval over which half the target initial mass is excavated 
by subcatastrophic collisions. Then, Me. = 0.5Mt,initial and 
when we integrate Eqn.[3]we obtain the timescale 

A . _ _ 12(s -I- 4)Mt, initial _ , 

where we used Ip = Dp7t/6 for spherical projectiles, s = 
—2.574 is the exponent of the distribution n{Dp) — BDp, Dp 
is the projectile diameter and we approximated (Dt-fDp)^ ~ 
Dt, since the projectile’s diameter is much smaller than the 
target’s diameter (few percent for the largest projectiles in 
our simulations). For s = —2.574 used in our simulations, 
the calculation of At is not very sensitive to specific val¬ 
ues of Dp,max and Dp,min, because we measure the erosion 
efficiency of every projectile that collided with the target 
asteroid. When we set Dp,min = 1 mm instead of 1 m used 
in our simulations, the calculated timescale barely changes. 
The constant A can be derived from the scaling laws or found 
from the graph of ejected mass vs. projectile diameter of all 
our simulations (see Fig. [T]). The constant B = 5.07km~°~^ 
for t he projectile size distribution we used (iBottke et al.l 
I 2 OO 5 II . More details on derivation of these relations are in 
Appendix El 

3.3 Number of collisions measure 

The second approach is based on a calculation of the number 
of collisions with specified population of projectiles for a 
single target. 

We fi rst calculate the colli sional lifetime of a model as¬ 
teroid as llFarinella et al.|[l998l l 

fshatter = 20 Myr^ft), (5) 

where Rt is the mean target asteroid radius in meters. We 
note that this collisional lifetime estimate is derived for 
smaller asteroids and it underestimates lifetime for aster¬ 
oids in our simulations. The collisional lifetime for aster- 
oids larger t h an se veral hundred meters as calculated by 
iBottke et al.l (l2005lj is more complicated function, but in 
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Table 3. Semiaxes ratios of the ellipsoid representing the target 
asteroid used in our simulations (semiaxes a, b, c). The shortest 
axis was c = 0.5 km. 


a/c 1.1 1.2 1.3 1.3 1.30 1.4 1.4 1.4 1.5 2.0 

b/c 1.1 1.2 1.3 1.2 1.14 1.4 1.1 1.2 1.2 1.4 



Figure 1. Erosion efficiency of projectiles in all our simulations 
measured by ejected mass vs. projectile diameter. The ejected 
mass is proportional to the volume of the projectile. The vertical 
spread of the points for a specific projectile diameter is caused by 
various impact geometries. 


genera l it has a steeper slope than that of iFarinella et ahl 
J 19981) . 

Then, we calculate a number of subcatastrophic col¬ 
lisions that the asteroid experiences during its collisional 
lifetime. In our simulations, we introduced a cutoff of the 
projectile size equal to 1 m because of quickly growing calcu¬ 
lation expenses with decreasing size of projectiles. We think 
that smaller projectiles are not very relevant for this effect 
but more extensive calculations will have to be done to check 
this in the future. 

The number of collis ions for a single target is 
llO’Brien fc Greenberdlioosh 

ricoii = Pi{Rt + Rp)^NpAt , (6) 


where pi is the intrinsic probability of collision, Rt and Rp 
are the target and the projectile radii, respectively, Np is 
the number of projectiles in this region and At is the time 
interval for which we calculate the number of collisions. We 
take this equal to the tshatter as calculated by Eqn. [5] We 
assume the model asteroid is a member of t he Main Asteroid 
Belt. T he intrinsic probability is taken from ICibulkova et al.l 
||2014| as Pi = 4.91 x 10"^® km”Vr"^- 

The number of projectiles is calculated assuming that 
the distribution of asteroid sizes in the Main Asteroid Belt 
is a power law as already mentioned above dBottke et al.l 
l2005h . The number of bodies in the size interval from Dp,min 

to Dp,max is 


V(Dp,min, Dp.max) = ^(D^+m',,, - , (7) 

where we use the same size distribution as in Eqns. [3]and[l 
and q = logjQS. The projectile sizes used for various target 


sizes are in T able [21 The upper limit is giv en by disrup - 
tion criterion dStewart fc Leinhardtl[201^ an dH ouse nl l2009lj 
and the lower limit is chosen somewhat arbitrarily to keep 
the computational expenses reasonable. Finally, a rough es¬ 
timate of a timescale of the increasing elongation can be 
calculated. It is 

, _ ^shatter ^reshape 

^reshape — 7 5 (oj 

'^coll Oavg 

where dreshape is the magnitude of a relative change of any 
semiaxis of the asteroid after that the asteroid will experi¬ 
ence some abrupt change of shape and 5avg is the average 
relative change of the same ratio per one collision. We dis¬ 
cuss the former value in the next paragraph, the latter value 
comes from our simulations. Also note, that this timescale 
actually does not depend on the asteroid lifetime (put the 

Eqn. [6] to the Eqn . [51). _ 

According to [Harris et al.l ([2009f) , there are no photo¬ 
metric observations of asteroids implying bodily axial ra¬ 
tio higher than about 3. Therefore, we chose the value of 
^reshape = 3/2 for the most elongated initial asteroid model 
used in our simulations (initial axial ratios 2 : 1.4 : 1) and 
^reshape — 2.7 for an almost spherical initial asteroid model 
(initial axial ratios 1.1:1.1:1). 


3.4 Comparison of the two measures 

We decided to prefer the timescale measure based on the 
ejected mass which is a kind of efficiency measure for every 
projectile. This measure gives more realistic timescale esti¬ 
mate and is also consistent with the relative timescale of the 
process based on our simulations. 

When we analyzed our simulations we concluded that 
substantial change of the axial ratio was only caused by the 
very largest projectiles (or the collisions close to the shat¬ 
tering threshold). Even though smaller projectiles are more 
populous thanks to the power-law distribution of their sizes 
(the exponent being lower than —2), their cumulative effect 
does not balance the effect of the largest projectiles. It is 
necessary to somehow weight the excavation efficiency of the 
projectiles otherwise the ricoii is overestimated and fshatter is 
consequently underestimated. Therefore the timescale mea¬ 
sure based solely on the number of all collisions is unrealis¬ 
tically skewed to the small values. 


4 RESULTS 

In our simulations, we observed that the elongation of the 
target asteroid gradually increases with increasing number 
of collisions. As described in Section [5] we did not calcu¬ 
late a change of the actual shape of the target but rather a 
change of the shape of a dynamically equivalent equal mass 
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Figure 2. Series of successive collisions onto two asteroids of 
various mean diameters represented by ellipsoids with initial 
a/c = 2.0 and fe/c = 1.4. The lower abscissa is cumulative pro¬ 
jectile mass divided by the initial target mass, the upper abscissa 
is cumulative excavated mass divided by the initial target mass. 
Changing axial rati os are in the o rdinate. We also plot the simple 
analytical model of iHarri 3 illOOOl l (black dots), the “weak 17.6- 
km” asteroid model with the same shattering threshold as the 
4.2-km asteroid model and we mark the collisional lifetime of the 
model asteroids. See text for details. 


ellipsoid. We calculated its semiaxes from the inertia ten¬ 
sor of the target asteroid gradually changed by consecutive 
collisions forming craters and hence removing material from 
the asteroid. 

An example of three series of collisions onto the tar¬ 
get asteroid with initial semiaxes ratios of a/c = 2.0 and 
b/c = 1.4 and mean diameters of 4.2 and 17.6 km, respec¬ 
tively, is shown in Fig. (2] The graph shows the evolution 
of a/c, and b/c axial ratios vs. cumulative projectile mass 
divided by the initial target mass on the lower abscissa and 
cumulative excavated mass divided by the initial target mass 
on the upper abscissa. It can be seen that the two ratios tend 
to increase as more projectiles impact the target. Other sim¬ 
ulations look similar to this example. 

In Fig. [2] we also plot three series of collisions for a 
model of a “weak 17.6-km” asteroid for which we set the 
same shattering threshold Qg as for the 4.2-km asteroid. 
We ran these series to evaluate our hypothesis that the ap¬ 
parently faster erosion of the larger asteroid models is an 
artifact caused by their higher shattering thresholds as will 
be discussed in Sect. [5] T he black dots in the graph de¬ 
pict the analytical model ofH arrisi (Il990h which represents 
the erosion by removal of the constant-depth layer from the 
ellipsoidal asteroid model. This model is size-independent 
(contrary to our model) and its apparent consistency with 
collisional erosion of our 17.6-km asteroid model is largely 
coincidental. 

In the following Fig.O there are results of the timescale 
calculation of the erosional process for various size of the 
target asteroid. It shows both, the time needed to erode 
half the initial target mass by subcatastrophic collisions off 
the asteroid (Sect. 13.21) and also the time of reshaping the 
target asteroid from the initial axial ratio a/c = 2.0 to the 
axial ratio a/c = 3.0 (Sect. [3)3ll . The timescale is in millions 
of years. It is obvious from Fig. |3]that both timescales are 


Figure 3. The timescale of the collisional erosion calculated in 
two ways. Orange circles are calculated as the time after which 
half the initial asteroid mass is ejected by collisions, dark-red 
squares are calculated using the number of small collisions expe¬ 
rienced by the target asteroid during its c ollisional lifetime. Th e 
collisional lifetime is calculated by Egn. l lFarinella et al.l|l998ll . 

much longer than collisional lifetime for most asteroids in 
our simulations. 

The displayed timescales calculated by two different 
methods are similar for smaller asteroids but they diverge 
for larger asteroids. As already discussed in Sect. 13.41 . we 
prefer the timescale calculation based on the ejected mass 
(orange circles in Fig. |3} and all our conclusions come from 
this calculation. 

The initial elongation of the asteroid (shorter series of 
collisions) leads to faster erosion. The timescale calculation 
is affected by larger mean diameter of an asteroid Dt in 
Eqn.|3]and also by wider interval of projectile sizes for more 
elongated asteroids. 


5 DISCUSSION 

The exp lanation for the results of our simulations can be 
found in iHarri 3 |I1990I) . Cratering erosion tends to change all 
dimensions by the same amount on average. Imagine a triax- 
ial ellipsoid subjected to isotropic random impacts. We will 
concentrate only on the changing dimensions of the longest 
and the shortest axes of the ellipsoid for a while. We observe 
that impacts erode relatively larger part of the shortest axis 
of the ellipsoid than its longest axis. This leads to the in¬ 
creasing elongation of the ellipsoid or the ratio of the longest 
to the shortest axes leghts. That is exactly what we see in 
our present simulations as shown in Fig. [S] 

However, from the results it is obvious that the 
timescale of this process is several orders of magnitude 
longer than the collisional lifetime of the asteroid. It is also 
visible in Fig.[3]that for increasing size of the target asteroid 
the erosion is faster. It is probably caused by the increasing 
specific energy threshold in gravity regime for the asteroid 
sizes from about 0.6 km used in our simulations (see, e.g., a 
review bv iHolsapple et ^l2002l and their Fig. 6). It means 
that larger asteroid can sustain larger collisions which exca¬ 
vate more material. In Fig [2] it can be seen that the “weak 
17.6-km” asteroid model behaves much like the 4.2-km aster- 
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Figure 4. Relative sizes of large craters observed on small solar 
system bodies vs. mean body radius. Relative size F is a ratio 
of the crater diameter and the mean body radius. Larger bodies 
usually bear relatively larger impact craters. Overlapping points 
are the same craters on the same objects. 


oid because it cannot sustain larger craters than this smaller 
asteroid. Therefore, we think that this strengthening (by 
gravity) of asteroids with increasing size is responsible for 
faster collisional erosion. 

This is also supported by direct observations of large 
impact craters on various small solar system bodies. We 
illustrate this in Fig. [4] which we comp iled from crate r 
sizes and the parent bodies sizes listed in Thom^ ( 1999 


lAsphauj ||200^ and [Burchell fc Leliwa-Kopvstvhskil (I201C 
The graph shows that larger bodies bear relatively larger 
impact craters. 

We came to the similar conclusion as KA03, however, 
our particular results differ substantially. Our models are 
quite different and we think it could be useful to combine 
their main features to study the role of small collisions in a 
more complex way. 

In our model, we did not work with the actual shape of 
the asteroid but with the dynamically equivalent ellipsoid 
as described above. Probably more serious drawback of our 
model is the absence of “landsliding” that was present in 
KAOS and is probably important for a shape relaxation to 
equilibrium (there is a limit of the angle of repose of the sur¬ 
face material observed on the asteroids studied by spacecraft 
mission as well as predicted by theory). 

On the other hand, our model deals with angular mo¬ 
mentum evolution of the target asteroid caused by projec¬ 
tiles, mass loss caused by impact excavation and we also al¬ 
low for large craters up to the shattering threshold. All these 
characteristics were mentioned by KAOS as future plans of 
the modelling of impact evolution of shapes of asteroids. 

Important difference between our and KAOS’s model is 
that we do not model ejecta’s orbits and their reimpact on 
the surface of the target asteroid. In our model, all ejected 
material escapes the asteroid. This makes the model only 
valid for small asteroids in the size range of several hundred 
meters to about 18-20km (see Sect. [2] for details). 

We also differ in the resulting shape of the asteroids 
subjected to a series of small collisions. While KAOS ob¬ 
tained oblate shape in most of their simulations, we pre¬ 
dict the elongation will grow for an already elongated initial 


shape. Our result seems to be consistent with an expected 
behaviour of large number of collisions as noted above. KAOS 
also noted the dependence of the result on the initial spin 
rate. We cannot verify this result since our model asteroids 
had random rotation period as described in Sect. [2] 

We predict even longer timescale of the erosion than 
KA03 which can be attributed to different construction of 
the model and the timescale calculation. However, the qual¬ 
itative result is the samethe timescale is several orders of 
magnitude longer than the collisional lifetime of the aster¬ 
oid for the projectile population exponent s = —2.574. 

What follows are few caveats concerning the construc¬ 
tion of our model. It has several points that are affected 
by large uncertainties. First, scaling laws for impact crater 
formation were developed for halfspace impacts, not for im¬ 
pacts onto curved surfaces and also we used earth rock ma¬ 
terial constants which may or may not be realistic. 

Second, the catastrophic disruption threshold comes 
from numerical simulations and scaling and is only one of 
many that can be found in literature as there is quite a large 
uncertainty in this. 

Third, there is a large uncertainty in population esti¬ 
mates of such small projectiles (and also targets), since we 
are well beyond present observation size limit. The power- 
law exponent is taken fro m the numerical mod el of the Main 
Asteroid Belt evolution of iBottke et al.l (l2005ll . Other values 
are possible and would lead to faster or slower evolution of 
the asteroid shape. 

Fourth, there is quite a questionable assumption about 
the isotropic distribution of projectiles. We assumed that 
the Main Belt asteroids have some moderate range in orbit 
inclination and also that the rotation axes of the target as¬ 
teroids are randomly oriented. This last assumption may be 
erroneous, though, becaus e small Main Belt astero ids may 
be trapped in Slivan state llVokrouhlickv et al.ll2003l l caused 
by the Yarkovsky-O’Keefe-Radzievskii-Paddack (YORP) 
effect which makes obliquity of their rotation axis gain spe¬ 
cific values. It is, however, not clear whether this mecha¬ 
nism also works for smaller, km-sized asteroids. The above 
mentioned paper indicates that asteroids smaller than about 
10 km may be able to leave the resonant spin state. 


6 CONCLUSIONS 

In a present paper we describe the effect of increasing elon¬ 
gation of a model asteroid represented by a homogeneous 
triaxial ellipsoid caused by consecutive small collisions. The 
changing shape of our model asteroid was indicated by 
changing intertia tensor and axial ratios of a dynamically 
equivalent equal mass ellipsoid. The effect can be explained 
by the fact that the shorter axis of the ellipsoid decreases 
relatively faster than the longer one for large number of 
isotropically distributed collisions. 

However, the estimated timescale of this process seems 
to be much longer than collisional lifetime of asteroids in 
the size range of several hundred meters to about 20 km. 
Therefore, this effect is probably not very important in for¬ 
mation of overall asteroid shapes and their evolution, unless 
some of our assumptions is unrealistic. For instance, if the 
exponent of the power-law projectile size distribution was 
smaller than s = —2.574, the erosion would be faster. It is 
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therefore important to set new observational constraints on 
the small projectile population so that we can obtain more 
accurate timescale estimates. 

There is a group of newly discovered asteroids called 
Barbarians after its first member, 234 Barbara, that are 
thought to be very old according to t heir polarimetric and 
spectral properties (|Tanga et al.l 120151 1. These observations 
indicate rather high content of calcium-aluminium-rich in¬ 
clusions l|Sunshine et al.ll2008h . Collisional erosion may pos¬ 
sibly play some role for larger long living asteroids, but our 
present model cannot deal with reacumulated ejecta effects, 
that are necessary to include for asteroids larger than about 
20 km in diameter. Further investigation of this area is nec¬ 
essary. 
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APPENDIX A: EJECTED MASS TIMESCALE 
DERIVATION 

Here we derive the relations for erosion timescale estimate 
based on the excavated mass efficiency of projectiles from 
the Sect. [321 Eroded mass is given by the collision prob¬ 
ability of a target asteroid with a diameter Dt and a pro¬ 
jectile with a diameter Dp, an efficiency of mass excavation 
by the projectile given by its volume Vp and a pr oportional¬ 
ity co nstant A (it is given by the scaling laws of iHolsapplei 
I 2 OO 3 II per infinitesimal time interval t times projectile size 
distribution n{Dp). The probability of collision is given by 
the cross-sectional area of the sum of Dt and Dp times the 
intrinsic probability of collisio n pi (basically the flux of pr o- 
jectiles per area per time, see lO’Brien fc Greenberg|l2005l ~) 

Me. = ipi(A -f DpfAVpn{Dp)Ppt. (Al) 

We plugin the projectile size distr ibution n{Dp) = BDp (we 
used that from iBottke et aD2005l '). a spherical projectile vol¬ 
ume Vp = DpU/G and we integrate over a time interval At, 
projectile size interval from Dp,min to Dp,max on the right- 
hand side of Eqn. \M\ We also approximate (Dt-fDp)^ « Dt 
as the projectile’s diameter is much smaller than target’s di¬ 
ameter (few percent for the largest projectiles in our simula¬ 
tions). We can arbitrarily choose the interval of integration 
of the left-hand side of Eqn. \M\ The reasonable choice is to 
put JMbt = 0.5Mt,initial because a catastrophic disruption 
of an asteroid is defined the same way. 

When we integrate Eqn. \M\ we can finally evaluate 
the time At after which the target asteroid subjected to 
consecutive collisions with a population of projectiles looses 
half its original mass as 
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At 


12 (s + 4) Aft,initial 


PinABDl{D 


S+4 

p,max 


- D 


s+4 \ 
p,min/ 


(A2) 


When we plugin Aft, initial ~ to Ean. lA2l we see, that the 
timescale is proportional to the target’s mean diameter Dt, 
inversely proportional to pi (smaller flux of projectiles leads 
to a longer timescale) and inversely proportional to A, or 
erosion efficiency of projectiles. 


This paper has been typeset from a lAtjjjX file prepared 
by the author. 



